Numerical study of dynamo action at low magnetic Prandtl numbers 
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We present a three-pronged numerical approach to the dynamo problem at low magnetic Prandtl 
numbers Pm. The difficulty of resolving a large range of scales is circumvented by combining Direct 
Numerical Simulations, a Lagrangian-averaged model, and Large-Eddy Simulations (LES). The 
flow is generated by the Taylor-Green forcing; it combines a well defined structure at large scales 
and turbulent fluctuations at small scales. Our main findings are: (i) dynamos are observed from 
Pm = 1 down to Pm = 10 -2 ; (ii) the critical magnetic Reynolds number increases sharply with 
P M as turbulence sets in and then saturates; (iii) in the linear growth phase, the most unstable 
magnetic modes move to small scales as Pm is decreased and a Kazantsev k 3 ^ 2 spectrum develops; 
then the dynamo grows at large scales and modifies the turbulent velocity fluctuations. 

PACS numbers: 47.27.eq,47.65.+a91.25w 
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stretching and folding of magnetic field lines by veloc- 
ity gradients overcome dissipation, dynamo action takes 
place above a critical magnetic Reynolds number R\i, 
with Rm = PmRv = UL/rj. Dynamo experiments en- 
gineering constrained helical flows of liquid sodium have 
been successful 0- However, these experimental setups 
do not allow for a complete investigation of the dynam- 
ical regime, and many groups have searched to imple- 
ment unconstrained dynamos 0. Two difficulties arise: 
first, turbulence now becomes fully developed with ve- 
locity fluctuations reaching up to 40% of the mean; sec- 
ond, it is difficult to engineer flows with helical small 
scales so that the net effect of turbulence is uncertain. 
Recent Direct Numerical Simulations (DNS) address the 
case of randomly forced, non-helical flows with magnetic 
Prandtl numbers from 1 to 0.1. Contradictory results 
arc obtained: it is shown in || that dynamo action can 
be inhibited for Pm < 1/4, while it is observed in Q 
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v . y v = — V7 5 + j x B + z/Y7 2 v + F (1) that the dynamo threshold increases as P M down to 

Pm ~ 0.3. Experiments made in von Karman geometries 
^2) (either spherical or cylindrical) have reached Rm values 
up to 60 Q. Also, MHD turbulence at low Pm has been 
studied in the idealized context of turbulent closures • 
In this context, turbulent dynamos are found, and the 
dependences of R C M upon three quantities are studied, 
namely Pm, the relative rate of helicity injection, and 
the forcing scale. An increase of ~ 20% in R C M is ob- 
served as Pm decreases from 1 to ~ 3 x 10 -5 . Recently, 
the Kazantsev-Kraichnan || model of (^-correlated ve- 
locity fluctuations has been used to study the effect of 
turbulence. It is shown that the threshold increases with 
the rugosity of the flow field Q , and that turbulence can 
either increase or decrease the dynamo threshold depend- 
ing on the fine structure of the velocity fluctuations [ljj • 



The generation of magnetic fields in celestial bodies 
occurs in media for which the viscosity v and the mag- 
netic diffusivity r\ are vastly different. For example, in 
the interstellar medium the magnetic Prandtl number 
Pm = v/r) has been estimated to be as large as 10 14 , 
whereas in stars such as the Sun and for planets such 
as the Earth, it can be very low (Pm < 10~ 5 , the value 
for the Earth's iron core). Similarly in liquid breeder 
reactors and in laboratory experiments in liquid met- 
als, Pm « 1. At the same time, the Reynolds number 
Ry = UL/v (U is the r.m.s. velocity, L is the integral 
scale of the flow) is very large, and the flow is highly com- 
plex and turbulent, with prevailing non-linear effects ren- 
dering the problem difficult to address. If in the smallest 
scales of astrophysical objects plasma effects may prevail, 
the large scales are adequately described by the equations 
of magnetohydrodynamics (MHD), 



v VB = B Vv + ?7V 2 B , 



together with V • v = 0, V • B = 0, and assuming a con- 
stant mass density. Here, v is the velocity field normal- 
ized to the r.m.s. fluid flow speed, and B the magnetic 
field converted to velocity units by means of an equiva- 
lent Alfven speed. V is the pressure and j = VxB the 
current density. F is a forcing term, responsible for the 
generation of the flow (buoyancy and Coriolis in planets, 
mechanical drive in experiments). 

Several mechanisms have been studied for dynamo ac- 
tion, both analytically and numerically, involving in par- 
ticular the role of helicity Q (i.e. the correlation between 
velocity and its curl, the vorticity) for dynamo growth 
at scales larger than that of the velocity, and the role 
of chaotic fields for small-scale growth of magnetic ex- 
citation (for a recent review, see [j|). Granted that the 



There is therefore a strong motivation to study how 
the dynamo threshold varies as Pm is progressively de- 
creased, for a given flow. In this letter we focus on a 
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situation where the flow forcing is not random, but gen- 
erates a well defined geometry at large scales, with tur- 
bulence developing naturally at small scales as the Ry 
increases. This situation complements recent numerical 
works 0, 0, EJ, E3 and is quite relevant for planetary and 
laboratory flows. Specifically, we consider the swirling 
flow resulting from the Taylor-Green forcing fll] |: 
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— cos(fco x) sin(fco y) cos(fco z) 
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with kg = 2, so that dynamo action is free to develop 
at scales larger or smaller than the forcing scale kf = 
ko^/3. This force generates flow cells that have locally 
differential rotation and helicity, two key ingredients for 
dynamo action 0, Note that the net helicity, i.e. 
averaged in time and space, is zero in the 27r-periodic 
domain. However strong local fluctuations of helicity are 
always present in the flow. Small scales are statistically 
non-helical. The resulting flow also shares similarities 
with the Maryland, Cadarache and Wisconsin sodium 
experiments and it has motivated several numerical 
studies at P M ~ 1 EES- 
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TABLE I: Parameters of the computation: code used, 
linear grid resolution N, Reynolds number Rv, integral 
length scale L (defined from the kinetic energy spec- 
trum L = 2-7T J fc _1 Ev(k)dk/ J Ev(k)dk), critical magnetic 
Reynolds number R C M , inverse magnetic Prandtl number 
1/P M , wavenumber /cmax with the largest magnetic energy, 
characteristic wavenumber ko of magnetic field gradients (de- 
fined as the maximum of the current density spectrum), and 
kinetic spectral index p in the range [fcMAx,M- The values 
of p, L and U used in the definitions of the Reynolds and mag- 
netic Prandtl numbers, are computed as time averages during 
the steady state of the hydrodynamic simulation; kMAx and 
ko are computed as time averages during the linear regime of 
the dynamo simulation closest to R M . 



Our numerical study begins with DNS in a 3D periodic 
domain. The code uses a pseudo-spectral algorithm, an 
explicit second order Runge-Kutta advance in time, and 
a classical dealiasing rule — the last resolved wavenum- 
ber is k = N/3 where N is the number of grid points per 
dimension. Resolutions from 64 3 to 512 3 grid points are 
used, to cover Pm from 1 to 1/5. However, DNS are lim- 
ited in the Reynolds numbers and the (lowest) Pm they 
can reach. We then use a second method, the LAMHD 
(or a) model, in which we integrate the Lagrangian- 
averaged MHD equations 0, [l^. This formulation 
leads to a drastic reduction in the degrees of freedom 
at small scales by the introduction of smoothing lengths 
ay and olm- The fields are written as the sum of filtered 
(smoothed) and fluctuating components: v = u s + 5v, 



B = B S + <5B, with u s = G a , 



B s = Gr 



)B, where 



'(2)' stands for convolution and G a is the smoothing kernel 
at scale a, G a (r,t) = exp[— r/a\/A-Ka 2 r. Inversely, the 
rough fields can be written in terms of their filtered coun- 



u s andB = (l-a 2 M V 2 ) B 



tcrparts as: v = (1- 
In the resulting equations, the velocity and magnetic field 
are smoothed, but not the fields' sources, i.e. the vortic- 
ity and the current density |l6j . This model has been 
checked in the fluid case against experiments and DNS 
of the Navier-Stokes equations in 3D E3> as well as in 
MHD in 2D ^3- Finally, in order to reach still lower 
Pm, we implement an LES model. LES are commonly 
used and well tested in fluid dynamics against laboratory 
experiments and DNS in a variety of flow configurations 



but their extension to MHD is still in its infancy 
(see however We use a scheme as introduced in |20j . 
aimed at integrating the primitive MHD equations with a 
turbulent velocity field all the way down to the magnetic 
diffusion with no modeling in the induction equation but 
with the help of a dynamical eddy viscosity [21| : 

v{k, t) = 0.27[1 + 3.58(fc/A' c ) 8 ] y/E v (K c ,t)/K c ; (4) 
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FIG. 1: R%j for dynamo action versus inverse Pm. Symbols 
are: x (DNS), + (LAMHD), and o (LES). Transverse lines 
indicate error bars in the determination of R M , as the distance 
between growing and decaying runs at a constant Rv. 
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FIG. 3: Spatial distributions of the magnetic energy for two 
Taylor-Green cells (DNS) : P M = 1, Rv ~ 30 at t = 20 (left), 
Pm = 0.4, R v = 675 at t = 20 (center), and t = 150 (right). 



FIG. 2: Magnetic spectra for P M = 1 to P M = 0.4 (DNS), 
Pm = 0.25,0.13 (LAMHD), P M = 0.038 (LES), at a time 
within the linear growth of magnetic energy. 

K c is the cut-off wavenumber of the velocity field, and 
Ev(k,t) is the one-dimensional kinetic energy spectrum. 
A consistency condition for our approach is that the mag- 
netic field fluctuations be fully resolved when 2ir/K c is 
smaller than the magnetic diffusive scale £ n ~ L/Rm- 

The numerical methods, parameters of the runs, and 
associated characteristic quantities are given in Table I. 
In all cases, we first perform a hydrodynamic run, lasting 
about 10 turnover times, to obtain a statistically steady 
flow. Then we add a seed magnetic field, and monitor 
the growth of the magnetic energy Em for a time that 
depends on the run resolution; it is of the order of 1 mag- 
netic diffusion time t v = (2tt) 2 /i] at 64 3 , but it drops 
down to t j; /5 at 512 3 . We define the magnetic energy 
growth rate as a = d log Eujdt, computed in the linear 
regime (t is in units of large scale turnover time). The 
dynamo threshold corresponds to a = 0. For each config- 
uration (Table I) , we make several MHD simulations with 
different Pm, varying r], and for a fixed Rv defined by 
the hydrodynamic run. We bound the marginal growth 
between clearly decaying and growing evolutions of the 
magnetic energy. This procedure is unavoidable because 
of the critical slowing down near threshold. 

At Pm = 1, the dynamo self-generates at R C M = 30. As 
Pm is lowered, we observe in the DNS that the thresh- 
old reaches R C M = 70 at Pm = 0.5 and then increases 
steeply to R C M = 226 at Pm = 0.3; at lower Pm it does 
not increase anymore, but drops slightly to a value of 
200 at P M = 0.2 (Fig.l and Table I). Wc then continue 
with LAMHD simulations to reach lower Pm- To en- 
sure the consistency of the method, we have run over- 
lapping DNS and LAMHD simulations in the range from 
Pm = 0.4 — 0.2, the agreement of the two methods being 
evaluated by the matching of the magnetic energy growth 
(or decay) rates for identical (Pm,Rm) parameters. We 
have observed that a good agreement between the two 
methods can be reached if one uses two different filtering 
scales av and a m in LAMHD, chosen to maintain a di- 
mensional relationship between the magnetic and kinetic 
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dissipation scales, namely ay /cum = Pm ■ ^ ur obser- 
vation with the LAMHD computations is that the steep 
increase in R° M to a value over 250 is being followed by a 
plateau for Pm values down to 0.09. We do note a small 
but systematic trend of the LAMHD simulations to over- 
estimate the threshold compared to DNS. We attribute 
it to the increased turbulent intermittency generated by 
the a model, but further investigations are required to 
describe fully this effect. The LES simulations allow us 
to further our investigation; with this model the thresh- 
old for dynamo self-generation remains constant, of the 
order of 150, for Pm between 10 _1 and 10 -2 . 

In regards to the generation of dynamo action in the 
Taylor-Green geometry we thus find: (i) at all Pm in- 
vestigated a dynamo threshold exists; (ii) as Pm drops 
below 0.2 - 0.3, the critical R C M levels and remains of 
the order of 200; (iii) the steep initial increase in R C M is 
identified with the development of an inertial range in the 
spectra of kinetic energy. As the kinetic energy spectrum 
grows progressively into a Kolmogorov fc -5 / 3 spectrum, 
R° M ceases to have significant changes - cf. Table H] 

We plot in Fig. |2] the magnetic energy spectra Em (k) 
during the linear growth phase, at identical instants when 
normalized by the growth rate. Four features are note- 
worthy: first, the dynamo grows from a broad range of 
modes; second, the maximum of Em{^) moves progres- 
sively to smaller scales as Pm decreases, a result already 
found numerically in Q; third, a self-similar magnetic 
spectrum, EM(k) ~ fc 3 / 2 , develops at the beginning dur- 
ing the linear growth phase — as predicted by Kazant- 
sev || and found in other numerical simulations of dy- 
namo generation by turbulent fluctuations 0, 01 • This is 
a feature that thus persists when the flow has well defined 
mean geometry in addition to turbulence. Lastly we ob- 
serve that the initial magnetic growth at small scales is 
always followed by a second phase where the magnetic 
field grows in the (large) scales of the Taylor-Green flow. 
Figure 3 shows renderings of the magnetic energy and 
compare low and high Reynolds number cases. When the 
dynamo is generated at low Reynolds number (Rv ~ 30 
and Pm = 1), the magnetic field is smooth. As Pm 
decreases and the dynamo grows from a turbulent field, 
one first observes a complex magnetic field pattern - for 



4 



M 10 




PCMI Programs for support. Computer time was pro- 
vided by NCAR, PSC, NERSC, and IDRIS (CNRS). 



FIG. 4: Magnetic (thin line) and kinetic (thick line) spectra 
as a function of time at Pm = 0.4 (DNS). 



t < 40, in the example shown in Fig. 3 (center). But as 
non-linear effects develop (here for times t > 40) a large 
scale mode (k = 2) dominates the growth with a struc- 
ture that is similar to the one at low Ry. The initial 
growth of small scale magnetic fields and the subsequent 
transfer to a large scale dynamo mode is also clearly vis- 
ible on the development in time of the magnetic and ki- 
netic energies, in a high Ry shown in Fig. 4. 
During the linear growth, a wide interval of modes in- 
crease in a self-similar fashion, accounting for the com- 
plexity of the dynamo field - cf. Fig. 3 (center). At a 
later time, the large scale field grows and the kinetic en- 
ergy spectrum Ey(k) is progressively modified at inertial 
scales. The spectral slope changes from a Kolmogorov 
A: -5 / 3 scaling to a steeper, close to fc -3 , regime (22}. The 
effect is to modify the turbulent scales and to favor the 
dynamo mode that is allowed by the large scale flow ge- 
ometry. This is consistent with the development of a 
k magnetic spectrum, observed in the Karlsruhe dy- 
namo experiment [2^] . ft also corroborates the claim j2J| 
that the saturation of the turbulent dynamo starts with 
the back-reaction of the Lorentz force on the turbulent 
fluctuations. 

To conclude, using a combination of DNS, LAMHD 
modeling and LES, we show that, for the Taylor-Green 
flow forcing, there is a strong increase in the critical mag- 
netic Reynolds number for dynamo action when Pm is 
decreased, directly linked to the development of turbu- 
lence; and it is followed by a plateau on a large range of 
Pm from ~ 0.25 to ~ 10~ 2 . In a situation with both a 
mean flow and turbulent fluctuations, we find that the 
selection of the dynamo mode results from a subtle in- 
teraction between the large and small scales. 
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